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We study the quantum dynamics of the cavity optomechanical system formed by a Fabry-Perot 
cavity with a thin vibrating membrane at its center. We determine in particular to what extent 
optical absorption by the membrane hinders reaching a quantum regime for the cavity-membrane 
system. We show that even though membrane absorption may significantly lower the cavity finesse 
and also heat the membrane, one can still simultaneously achieve ground state cooling of a vibra- 
tional mode of the membrane and stationary optomechanical entanglement with state-of-the-art 
apparatuses. 
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I. INTRODUCTION 



Cavity optomechanics has recently emerged as a 
rapidly developing research held in which relevant tasks 
such as quantum limited displacement sensing, highly 
sensitive force measurements, and the realization of 
quantum interfaces for quantum information networks, 
may be realized by virtue of the radiation pressure cou- 
pling between the light confined in a cavity and nano- and 
micro- mechanical resonators [IHl]- The standard and 
simplest optomechanical setup is a Fabry-Perot cavity 
in which one of the two mirrors is a vibrating microme- 
chanical object. It has been the first setup experimen- 
tally studied 04^]: and for which ground state cooling 
of the resonator has been recently approached Q, and 
strong optomechanical coupling has been demonstrated 
[Toj | . In order to enter in a regime where the dynamics 
is distinctly quantum, the optomechanical coupling has 
to dominate both optical and mechanical losses, which 
implies a high-finesse optical cavity enabling to amplify 
the radiation pressure, and a movable mirror with high 
mechanical quality factor Q m . However, it is technically 
non-trivial to fabricate devices possessing simultaneously 
high Q m and high reflectivity, and therefore novel devices 
have been recently designed and tested. A first exam- 
ple is provided by silica toroidal optical microresonators, 
coupled by the radiation pressure of the circulating light 
with the radial vibrational modes of the support [l|, llll - 
lill . Other novel approaches propose to avoid completely 
clamping mechanical losses by optically trapping levitat- 
ing dielectrics as in [l4j]. Furthermore, other optome- 
chanical devices have been tested, in which the coupling 
is provided by the transverse gradient force, such as sus- 
pended silicon photonic waveguides [l5[, SiN nanowire 
evanescently coupled to a microtoroidal resonator 
and a "zipper" cavity formed by two adjacent photonic 
crystal wires (l7j . 

A further solution has been proposed in Refs. [Ill [l9[ 
in which the mechanical degree of freedom is represented 



by a partially reflective thin vibrating membrane inserted 
within a standard cavity with two bulk mirrors. A high 
value of Q m is reached by using commercially available 
high-stress SiN membranes, which are known to possess 
very low intrinsic losses [20(. At the same time one has 
a high-finesse cavity by using standard highly reflective 
mirrors and minimizing optical absorption in the mem- 
brane, which can be achieved using stoichiometric SiN 
pH ] and very thin membranes. As recently shown in 
Ref. [IH, this "membrane-in-the-middle" setup has the 
further advantage that the radiation pressure coupling 
can be tuned in situ by adjusting the position and orien- 
tation of the membrane within the cavity [22j . In partic- 
ular optomechanical couplings which are nonlinear in the 
resonator position can be engineered and tuned. In this 
paper we study this cavity-membrane system by investi- 
gating to which extent optical absorption by the mem- 
brane affects reaching the quantum regime. We shall sec 
that even if absorption appreciably decreases the cavity 
finesse and also heats the membrane, it is still possible to 
find realistic parameter settings where quantum behavior 
of the cavity-membrane system is detectable, even up to 
room temperature. 



In Sec. II we describe the system by means of quan- 
tum Langevin equations, showing in particular how the 
dynamics is modified by optical absorption by the mem- 
brane. In Sec. Ill we focus on the quantum behavior 
of the optomechanical system, and show that ground 
state cooling of the mechanical resonator, and steady- 
state optomechanical entanglement are achievable even 
for a non-negligible absorption. Sec. IV is for conclud- 
ing remarks, while details of the derivation of the stan- 
dard single-mode cavity optomechanics description from 
the general multimode radiation pressure interaction are 
given in the Appendix. 
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II. DESCRIPTION OF THE SYSTEM 

We consider a thin semi-transparent dielectric mem- 
brane placed within a driven high-finesse Fabry-Perot 
cavity of length L (see Fig. [I]). The membrane vibra- 
tional modes are coupled to the optical cavity modes by 
radiation pressure, and therefore one has in general a 
multimodc bosonic system in which mechanical and op- 
tical modes interact in a nonlinear way. However, in most 
practical situations, one can adopt a simplified descrip- 
tion based on a single cavity mode interacting with a 
single mechanical mode [l8|, [l9j| . One can restrict to a 
single cavity mode whenever the driving laser populates 
a given cavity mode only (typically a TEMoo mode, here 
associated with the annihilation operator a) , and if scat- 
tering from the driven mode to other modes is negligible. 
This latter condition is satisfied when the frequencies of 
the membrane vibrational modes fli are smaller than the 
typical cavity mode separation, which is of the order of 
the free spectral range c/2L (see also Ref. [HI). More- 
over, one can consider a single mechanical mode of the 
membrane (described by dimensionless position q and 
momentum p operators, such that [q,p] = i) when the 
detection bandwidth is chosen so that it includes only 
a single, isolated, mechanical resonance with frequency 
Cl m . In Appendix |21 the general multimode description 
of the membrane-in-the-middle setup is provided and we 
discuss the conditions under which such a single-mode 
reduction is possible. 



included the radiation pressure interaction within the 
cavity energy term by defining the position-dependent 
frequency 

uj(q) =uj + 5oj[z Q (q)}, (2) 

where loq is the cavity mode frequency in the absence 
of the membrane, and 5ui[zo(q)] is the frequency shift 
caused by the insertion of the membrane (see Eq. (|A9|) 
for its explicit form and Appendix A for its derivation). 
This shift depends upon the membrane position along 
the cavity axis zo(q), which in turn depends upon the 
coordinate q because one can write zo(q) = z$ + Xaq, 
where zq is the membrane center-of-mass position along 
the cavity axis (see Fig. 1) and xq = \/h/mflm, with m 
the effective mass of the mechanical mode [l8|, [l9|, [2l[ ■ 

A first order expansion in q of Eq. © provides an accu- 
rate description of the physics in most cases, except when 
the membrane center zq is placed exactly at a node or at 
an antinode of the cavity field, where the first-order term 
in the expansion of oj(q) vanishes, and one has to consider 
the higher-order term, which is quadratic in q. This latter 
term describes a dispersive interaction between the opti- 
cal and the vibrational modes which has been discussed 
in [III [22|, HH [25[ and may be exploited for a quantum 
non-demolition measurement of the vibrational energy. 



A. Quantum Langevin description of dissipation 
and noise sources 



driving laser 



z = ok 



membrane 



The dynamics are determined not only by the Hamil- 
tonian of Eq. ([T]), but also by the fluctuation-dissipation 
processes affecting both the optical and the mechanical 
mode. The mechanical mode is affected by a viscous force 
with damping rate 7 m and by a Brownian stochastic force 
with zero mean value £ (t) , obeying the correlation func- 
tion siii 



L/2 



FIG. 1: (Color online) Schematic description of the cavity 
formed by two identical mirrors, with a thin semi-transparent 
vibrating membrane placed inside. 
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(3) 

where ks is the Boltzmann constant and Tq is the tem- 
perature of the reservoir of the membrane. The Brow- 
nian noise £(t) is a non-Markovian, Gaussian quantum 
stochastic process, but in the limit of high mechani- 
cal quality factor Q m = Q m /~/ m 3> 1, becomes with a 
good approximation Markovian, with correlation func- 
tion Salll 



By explicitly including cavity driving by a laser with 
frequency cjj, and input power V, one ends up with the 
following cavity optomechanical Hamiltonian, 



H = 



■ hu)(q)aJa + ihE(a 



(1) 

where E = y^T^Kn/fi^L, with kq the cavity mode band- 
width in the absence of the membrane. In Eq. (TTJ) we have 



{2n + l)5(t - t') + i 



S'(t-t') 



, (4) 



where no — (exp{ftr2 m /fcsTo} — 1) is the mean ther- 
mal phonon number at To, and S'(t — t') denotes the 
derivative of the Dirac delta. 

Then, due to the nonzero transmission of the cavity 
mirrors, the cavity field decays at rate kq and is affected 
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by the vacuum optical input noise a™ (tV whose only 
nonzero correlation functions is given by |27| 



(5) 



However, the presence of the membrane and of its op- 
tical absorption provides an additional loss channel for 
the cavity photons, which has been neglected up to now 
[H|,[24[. Absorption by the membrane is quantified by the 
imaginary part to the the refraction index, Hm = TiR+ini 
(ni <C nil), yielding a small imaginary part of the fre- 
quency shift Su> [zo(q)]. As a consequence, the membrane 
absorbs cavity photons with a rate 



Ki[q) = \lm{Suj [zo(q)]}\ , 



(6) 



which depends upon the effective resonator position op- 
erator q and therefore describes a nonlinear dissipative 
process affecting also the mechanical mode. Due to 
fluctuation-dissipation theorem, this nonlinear dissipa- 
tion is associated with an additional optical input noise 
a^ n (t), possessing the same autocorrelation functions of 
dg n (t) [see Eq. ((5J], and describing the vacuum optical 
field fluctuations re-radiated within the cavity by the 
membrane. At the same time this latter noise is respon- 
sible for an additional stochastic force on the mechanical 
mode, describing membrane heating due to absorption. 

Therefore, the cavity mode has a total decay rate 
k t(<?) = + K-i(q); this means that the cavity finesse 
is decreased from its empty cavity value J- by membrane 
absorption, and that, due to Eq. ([5]), it becomes depen- 
dent upon the membrane position Zq according to 



1 



1 



T T (z ) T 



Im 



LSuj(zq) 



(7) 



Using the explicit expression of the frequency shift de- 
rived in Appendix A, Eqs. (|A9[) . Eq. ([7]) implies that the 
resulting finesse Ft{zq) is a periodic function of zq with 
period equal to half wavelength of the mode, with max- 
ima attained when the membrane is centered around a 
node of the cavity field [H, [2l[. Such a behavior, in 
the case of an empty cavity finesse T = 20000 and for 
ni = 10 -4 , is shown in Fig. [2 where Ft{zq) decreases 
down to almost 50% of its initial value by varying the 
membrane position zq. 

Adding the above damping and noise terms to the 
Heisenbcrg equations of motion associated with the 
Hamiltonian of Eq. (TTJ) , one gets the following set of non- 
linear quantum Langevin equations (QLE) which, in the 
frame rotating at the driving laser frequency ujl, read 



q = f2 TO p, 

p = -Q m q - j m p - d q ui{q)a) a 



dqKxiq) 



V2ki(?) 

i [uj(q) - u>l\ a - [n + Ki(q)] a + E 



(8a) 
(8b) 

(8c) 
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FIG. 2: Cavity finesse versus the membrane position along 
the cavity axis zo- The plot refers to an empty cavity finesse 
J- = 20000, and to a membrane with thickness Ld = 50 nm, 
and index of refraction um = 2.0 + i!0 -4 . 



where d q denotes the derivative with respect to q. These 
nonlinear QLE exactly describe the dynamics of the 
cavity-membrane system, but they are difficult to solve 
due to the various nonlinearities. Membrane absorption 
in particular is responsible for the additional nonlinear 
damping term acting on the cavity mode, and for the ad- 
ditional nonlinear noise term associated with membrane 
heating. 



B. Linearized quantum Langevin equations 

In most practical situations it is not necessary to solve 
the nonlinear QLE of Eqs. ([S]). In fact, one is inter- 
ested in parameter regimes where the radiation pressure 
interaction has appreciable effects and a simple way to 
achieve it is to have an intense intracavity field, which is 
achievable with a large cavity finesse T and with enough 
driving power. In this case the operation point is given 
by a classical steady state characterized by a coherent 
intracavity field with amplitude a s (\a s \ 3> 1) and a de- 
formed membrane mode with a new stationary position 
q s , satisfying the coupled nonlinear conditions 



d q uj(q s )\a s 



\E\ 



[k + ki(<? s )] +[u} L -u(q s )Y 



(9a) 
(9b) 



Eq. (|9b[) in particular is the nonlinear equation determin- 
ing |a s |, due to the dependence of q s upon |a s | itself given 
by Eq. (|9aj) . and which may show optical bistability, i.e., 
the presence of two simultaneous stable solutions above a 
given threshold for the input power V, as demonstrated 
in [29[ and also with Bose-Einstein condensates in (30j . 

In this regime, the relevant dynamics concern the 
quantum fluctuations of the cavity and mechanical 
modes around the classical steady-state described above. 
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Rewriting each Heisenberg operator of Eqs. (O as the 
classical steady state value plus an additional fluctuation 
operator with zero mean value, and neglecting all the 
nonlinear terms in the equations, one gets the following 
linearized QLE for the fluctuations 



Sq 
Sp 



SX 



SY 



Sl m 8p, (10a) 
- [n m + d 2 q u{q s )\a s \ 2 } Sq - lm Sp + GSX 



(10b) 

V KIWI,) 

k-o + Ki(q s )] SX + A5Y - V2d q Ki(q s )a s 5q 

(10c) 



+V2^^o" + v^ites)xr, 

- [«o + Kx{q s )} SY - ASX + GSq 
+V^Y Q m + y/2 Kl (q a )Yr. 



(lOd) 



These equations have been derived by choosing the phase 
reference of the cavity field so that a s is real and posi- 
tive, and by defining the detuning A = uj(q s ) — ^l- Wc 
have also introduced the cavity field amplitude quadra- 
ture SX = (6 a + 8a [ )/\/2, the phase quadrature SY = 
(8a — 8a))/iy/2 and the corresponding Hermitian input 
noise operators Xf 1 = (of 1 + aJ lA )/y/2 and Yj n = 
(af - af A )/iV2, j = 0,1. The linearized QLE of 
Eqs. (|10[) show that the mechanical and cavity mode fluc- 
tuations are coupled by the effective optomechanical cou- 
pling 



G = -d q u(q s )a s V2 
= -2d q uj(q s ) 



(11) 



hbJr 



[K Q + Kl (q s )} 2 + A< 



which can be made very large by increasing the intracav- 
ity amplitude a s . 

The linearized QLE of Eqs. (fTUf provide a complete 
description of the dynamical effects of optical absorption 
by the membrane. These QLE are identical to the set 
of QLE describing standard cavity optomechanical sys- 
tems (see e.g. Refs. @, [U |32[) except for some addi- 
tional damping and noise terms: i) the cavity decay and 
noise terms related to Ki(q s ), describing the photon ab- 
sorbed by the membrane and the vacuum optical field re- 
emitted within the cavity; ii) the noise term proportional 
to Y™ describing heating of the mechanical resonator; iii) 
the dissipative coupling term y/2d q Ki(q s )a s 8q. Actually, 
in Eq. (|10b[) there is a further nonstandard term, given 
by — d q oj(q s )\a s \ 2 5q, describing a shift of the mechanical 
frequency caused by the nonlinear dependence of ui(q) 
upon q, which is not related with absorption. This latter 
term should be considered even when the quadratic term 
in oj(q) is very small, because it is proportional to the 
stationary number of intracavity photons |a; s | 2 , which is 
necessarily very large in the linearized regime we are con- 
sidering. 



III. EFFECT OF MEMBRANE ABSORPTION 
ON GROUND STATE COOLING AND 
STATIONARY ENTANGLEMENT 

As discussed in detail in the literature (3, HTT -ili . the 
radiation pressure interaction may lead a cavity optome- 
chanical system into a stationary state with distinct 
quantum features, and specifically with the mechanical 
resonator very close to its quantum ground state, and 
showing robust continuous variable (CV) entanglement 
between the cavity and the vibrational modes. It is there- 
fore interesting to verify if and when absorption by the 
membrane may hinder achieving the quantum regime for 
the cavity-membrane system under investigation. 

Eqs. (|10[) can be rewritten in compact form as u(t) = 
Au(t)+n(t), where u T (t) = [Sq(t),Sp(t),5X(t),SY(t)\ 
is 



the vector of fluctuation operators 



(*) 

y^MYl n (t)] 



the noise vector, and the matrix A is the drift matrix 
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(12) 



where we have defined T 
d 2 Mls)\a s \ 2 . 



V2d q Ki(q s )a s , and h = 
. We are interested in the stationary state 
of the system, which is reached for t — > oo for ev- 
ery initial state if the system is stable (see Appendix 
B for the stability conditions, which we will consider 
to be always satisfied from now on). Due to the lin- 
earization and to the Gaussian nature of the noise op- 
erators, the steady state is a zero-mean Gaussian state, 
completely characterized by its covariance matrix (CM). 
The latter is given by the 4x4 matrix V with ele- 
ments Vi m — ( u i (oo) u m (oo) + u rn (oo) ui (oo)) /2, where 
u m (oo) is the asymptotic value of the m-th component 
of the vector u(t). Using standard methods [HIHII], one 
gets that the steady state CM V is given by the solution 
of the Lyapunov equation 



AV + VA T = -D, 



(13) 



with A given by Eq. ([12]) and where D is the 4 x 
4 diffusion matrix, which is related with the cor- 
relation functions of the noise vector elements by 
((n k (s)ni(s') + ni(s')n k (s))) /2 = D ki S(s - s'), and 
whose explicit expression is 








7m (2n 





1) 


r 

2 



4«i(g») 






Kr(q s ) 




(14) 



nrtes) J 



Equation (fTB"]) is a linear equation for the CM V which 
can be solved in a straightforward way, but whose solu- 
tion is very cumbersome and will not be reported here. 
The CM contains all the information about the steady 
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state: in particular, the mean energy of the mechanical 
resonator is given by 



U 



m 2 ) + m\ 



[Vn + V 22 ] (15) 



h£l m I n + - 



where n is the effective mean vibrational number of the 
selected membrane mode. Obviously, in the absence of 
radiation pressure coupling it is n = no, where no cor- 
responds to the actual temperature of the environment 
To. The optomechanical coupling with the cavity mode 
can be used to engineer an effective bath of much lower 
temperature T <C To, so that the vibrational mode is 
cooled. Moreover from the covariance matrix V one can 
also calculate the entanglement between the cavity and 
the vibrational modes and at the steady state. We adopt 
as entanglement measure the logarithmic negativity En, 
which is a convenient entanglement measure, easy to 
compute and also additive. It is defined as [37j 



E N = max[0, — ln2n ], 



(16) 



with 



V 



2 -i/2 



1 /2 

£00 - [£(V) 2 -4dety] 1/2 l (17) S abs {oj) 



S(F) ee det Vi + det V 2 - 2 det V c , 



(18) 



where V\ , V 2 and V c are the 2x2 sub-block matrices of 
V 



V 



Vi v c 

2 



(19) 



Determining V from the Lyapunov equation of Eq. (|13|) 
and using Eqs. (|15[) - (fT8"|) , one can therefore see if mem- 
brane absorption affects reaching a quantum regime for 
the cavity-membrane system. 

At first sight, such absorption may be very harmful to 
quantum behavior, due to the additional damping and 
noise terms acting on both the optical and mechanical 
mode. In order to have an intuitive picture of how pho- 
ton absorption by the membrane manifests its effects 
we focus on ground state cooling, i.e., on the evalua- 
tion of the steady state elements of the CM, Vn and V 22 



[see Eq. ([IS])]. Following the approach of Refs. [H Hj|, 
one can equivalently determine the stationary CM, V, 
by solving the linearized QLE in the frequency domain 
rather than in the time domain. For the position and mo- 
mentum variances of the mechanical resonator one gets 
(see Ref. [H) 

<«■>-££«■<">■ <*> 

(20) 

where the position spectrum can be written as 

S q (uj) = IXcffMI 2 [S t h(u) + S rp (u) + S a bs(u)}, (21) 

with 



Sth(w) = -tt— coth 1 



fl m \2k B T 
the thermal noise spectrum, 

G 2 K T {q s ) [A 2 + K 2 T {q s ) + UJ 2 



7m (2n + 1) (22) 



S rp (uj) 



[K 2 T {q s ) + ( W - A)2] [„2 (g s ) + (OJ + A) 2 } 



(23) 



the radiation pressure noise spectrum, and finally 

T 2 TGA [A 2 + K 2 T (q s ) - lo 2 ] 



4 Kl (q s ) [K 2 T (q s ) + (w - A) 2 ] [K 2 T (q s ) + (u> + A) 2 } 



(24) 

the additional noise spectrum associated with membrane 
absorption, containing both the term quadratic in T de- 
scribing membrane heating, and the term linear in T, as- 
sociated with the nonzero correlation between such heat- 
ing noise and cavity phase noise [see also the diffusion 
matrix of Eq. (fT4")) ]. 

A further effect of optical absorption by the membrane 
is the modification of the effective mechanical suscepti- 
bility XeS (w) appearing in Eq. (|2~Tj) . which is given by 



Xeff(w) 



n 2 -uj 2 -tLWY Gn m [GA-r[KT(g.)-*4 



(25) 

where we have defined = VL 2 m + hfl m . 

The latter can be read as the susceptibility of an oscil- 
lator with effective (frequency-dependent) resonance fre- 
quency and damping rate respectively given by 



GQ m {GA [4(g s ) 



A 2 ] ~TK T (q s ) [ K 2 T (q s ) + UJ 2 + A 2 ]}' 



7™(w) = 7m + 



+ (w - A)2] [4( 9s ) + ( w + A)2] 
Gn m {2GAn T {q s ) - T [K 2 T (q s ) + lu 2 - A 2 ] } 



[4 (q s ) + (w 



A)2][4(q s ) + ( w + A)2 



(26) 
(27) 



As discussed in [39(, in the limiting case when the me- chanical quality factor Q m is large enough, and the op- 
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tomechanical coupling is not too strong, G <C Sl m (which 
is the relevant regime for ground state cooling) , this effec- 
tive susceptibility XeS (w) remains peaked at f2 m , and one 
can approximately calculate the variances of Eqs. (|2T)|) by 
evaluating the noise spectra just at the peak frequency 
f2 m . The integral are then easily performed and one gets 



1 



1 7m n + + A+ 



with 



_ G 2 k t + GT(A ± n m ) 

± 4( 9fl ) + (A±n m )2 ■ 



(28) 
(29) 



It is straightforward to verify that A± are just the scat- 
tering rates of laser photons into the antiStokes (A—) 
and Stokes (j4_|_) sidebands centered at wl ^ fi m , re- 
spectively. Therefore A^ — A + is the net laser cooling 
rate and it is related to the effective mechanical damp- 
ing by the relation 7^f (Q m ) = 7,„ + A_ - A + |!38j. The 
scattering rates of Eq. (|2"9")l are the sum of two terms, 
corresponding to the scattering of photons through the 
cavity (the first term) or through membrane absorption 
(the second term). Eq. (f2"9"| generalizes the expression 
of these rates obtained in Refs. [H, 0, HI] to the case 
with non-negligible absorption. However the relevant re- 
sult of the present analysis is that of Eq. (|2"5)l , clearly 
showing that one can achieve ground state cooling when 
A_ A+ , 7 m and both no and the term associated with 
membrane heating, T 2 /8Ki(q s )j m , are not too large. 

A simple description of the behavior of the logarithmic 
negativity is more difficult to find, but nonetheless one 
may reasonably expect that when the above conditions 
for achieving ground state cooling are met, one could also 
achieve stationary optomechanical entanglement as well. 
More generally speaking, the above analysis of ground 
state cooling suggests that it is still possible to find a pa- 
rameter regime where quantum behavior of the station- 
ary state of the membrane-in-the-middle setup is achiev- 
able, even in the presence of non-negligible absorption. 

We have verified these predictions by considering pa- 
rameters similar to the experiment of Ref. [2l[ . We have 
verified in particular that the approximate expression 
of Eq. (|28[) satisfactorily describes the vibrational oc- 
cupancy n at not too large input power. To be more 
specific we have taken a SiN membrane with refraction 
index n m = 2.0+zlO -5 , side length 0.5 mm and thickness 
Ld = 50 nm, implying an effective mass m = 8.5 ng; we 
also assume a mechanical resonance frequency fi m /27r = 
10 MHz, a mechanical quality factor Q m = 6 x 10 6 , and a 
cavity with length L = 0.74 mm and finesse J- = 23000. 
The results are shown in Figs. [SKI where we plot si- 
multaneously the mean vibrational occupancy n and the 
logarithmic negativity En as a function of various pa- 
rameters. All figures show that despite absorption may 
appreciably affect the cavity finesse, one can choose an 
operation point where the cavity-membrane system can 
be put in a steady state with distinct quantum features, 



i.e., with a vibrational mode close to its ground state and 
a significant optomechanical entanglement. 

In Fig. [3] n and En are plotted versus the cavity de- 
tuning A (at fixed reservoir temperature To = 1 K and 
laser input power V = 28.5 mW) and one can see that, 
as known in the literature [3l|-|34j, one has significative 
simultaneous cooling and entanglement when A ~ f2 m , 
i.e., when the cavity is resonant with the anti- Stokes side- 
band of the driving laser. Fig. 0] shows n and En versus 
temperature To (at fixed detuning A = f2 m input power 
V = 28.5 mW) and, as expected, they both worsen for 
increasing temperature; it is remarkable however that 
quantum features n < 1 and En > persists up to 
room temperatures, To ~ 300 K in this optimal parame- 
ter choice. In Fig. [5] instead we keep A = Vl m and To = 1 
K fixed and we plot n and En versus the input power 
V . Increasing V means increasing the effective optome- 
chanical coupling G and therefore also optomechanical 
entanglement increases; on the contrary n is minimum 
at not too large input power and therefore at not too 
large coupling G, because for larger G higher-order scat- 
tering process tend to heat the vibrational motion (see 
Refs. [32111 Hl|). Finally in Fig. ©we study the depen- 
dence of n and En upon the membrane thickness Ld, by 
keeping all the other parameter fixed (A = f2 m , To = 1 
K, V = 28.5 mW). We notice that there are two interest- 
ing parameter regions showing good simultaneous cool- 
ing and entanglement: one at small thickness Ld ~ 30 
nm, and a second one at larger thickness, Ld — 150 nm. 
Small membrane thickness is obviously advantageous for 
quantum effects because it minimizes both mass m and 
optical absorption effects; however one has appreciable 
cooling and entanglement also for a thicker membrane 
because in this case, even if the mass m and absorp- 
tion are larger, one has maximum membrane reflectivity 
TZ(k ). In fact Eq. £52} shows that TZ{k ) is maximum 
when riMkoLd = 7f/2, i.e. for Ld = A/4njw, which corre- 
sponds to Ld — 130 nm for a laser wavelength A = 1064 



IV. CONCLUSIONS 

We have analyzed the cavity optomechanical sys- 
tem formed by a Fabry-Perot cavity with a thin semi- 
transparent membrane placed at its center. In particular 
we have investigated the effect of optical absorption by 
the membrane, which is usually neglected in standard 
treatments. Membrane absorption has an appreciable 
effect on the cavity finesse and we have studied if and 
when it may also hinder achieving a quantum stationary 
state of the cavity-membrane system. Such absorption is 
responsible for a modification of the scattering rates of 
laser photons into the Stokes and antiStokes sidebands, 
and for an effective increase of the initial membrane tem- 
perature due to absorption-induced heating. Considering 
state-of-the-art experimental parameters, we have seen 
that one can always optimize the operation point so that 
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FIG. 3: (Color online) Effective mean vibrational number n 
(blue, full line) and logarithmic negativity En (red, dashed 
line) versus the scaled cavity detuning A/r2 m , at fixed in- 
put power V — 28.5 mW and temperature To = 1 K. Other 
parameters are: cavity length L — 0.74 mm, cavity finesse 
J- = 23000, membrane thickness Ld = 50 nm, membrane 
side-length 0.5 mm, effective mode mass m = 8.5 ng, mechan- 
ical resonance frequency fi m /27r = 10 MHz, and mechanical 
quality factor Q m = 6 x 10 6 . 
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FIG. 4: (Color online) Effective mean vibrational number n 
(blue, full line) and logarithmic negativity En (red, dashed 
line) versus temperature at fixed input power V = 28.5 mW 
and detuning A = Q m . Remarkably, quantum features persist 
up to room temperature. The other parameters are as in 
Fig. El 



membrane absorption does not pose serious limitations 
for reaching quantum effects. Remarkably, using param- 
eter analogous to those of Ref. [2l|, ground state cool- 
ing and stationary optmechanical entanglement may per- 
sist up to room temperatures. The fact that stationary 
optomechanical entanglement is possible even with non- 
negligible absorption is particularly relevant: appreciable 
optomechanical entanglement is achieved only at large 
enough intracavity power: in this regime however mem- 
brane heating by absorption becomes relevant and it is 
not obvious that optomechanical entanglement is robust 
even against such heating. Stationary entanglement is 




FIG. 5: (Color online) Effective mean vibrational number n 
(blue, full line) and logarithmic negativity En (red, dashed 
line) versus input power at fixed temperature To = 1 K and 
detuning A = Q m . The other parameters are as in Fig. [3] 
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FIG. 6: (Color online) Effective mean vibrational number n 
(blue, full line) and logarithmic negativity En (red, dashed 
line) versus membrane thickness at fixed temperature To = 1 
K, input power V = 28.5 mW, and detuning A = Q m . The 
other parameters are as in Fig. [3] 



particularly important because it lasts for ever and can 
be extremely useful for storing and redistributing quan- 
tum information, or realizing quantum interfaces [40l. l4lj . 
Finally, the membrane-in-the-middle setup could be ex- 
ploited also for quantum limited detection of weak forces 

Note added in proof Recently Ref. [43| derived the 
multimodc Hamiltonian of the one-dimensional limit of 
the membrane- in-the- middle scheme, without making the 
linear approximation in the interaction. The multimode 
Hamiltonian derived here in Appendix A coincides with 
that of Ref. [43| when the one-dimensional and linearized 
limit is taken. 
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[44l | . The cavity waist wo is related to the cavity geom- 
etry by the relation w 2 = {L/k)^J{l + g)/(l — g), where 
g = 1 — L/R, and the Rayleigh range zr is given by 
zr = kw 2 /2. Each mode oscillates in time with an angu- 
lar frequency w^ np given by 



Appendix A: General multi-mode Hamiltonian of 
the system 

An empty Fabry-Perot cavity supports an infinite set 
of optical modes, well described by the Hermite-Gauss 
modes [44| ■ In the case of a symmetric cavity formed by 
two identical spherical mirrors with radius of curvature 
R, separated by a distance L, in the coordinate system 
with the z axis along the cavity axis and centered at the 
cavity center (see Fig. 1), these Hermite-Gauss modes are 
written as 



^mnpt^i V: — T-mnpi^ Hi %) §m $ mn p(^, 7/, z) 



pn 



where, for to, n = 0, 1 . . ., p = 1, 2, . . ., 



(Al) 



w(z) 



w(z) 



exp 



x_+y_ 



w(z)y / n2 n + m - 2 nlm\L 



(A2) 



x ~t~ y 

®mn P {x, V, z) =kz- lp(z)(m + 71 + 1) + k ' 2 R( Z ) ' 

Here H n {x) is the n-th Hermite polynomial, w(z) = 
+ z 2 / z^ is the width of the cavity mode at z, 
R(z) = z + z 2 R j 1 z is the radius of curvature of the wave- 
front at z, ip( z ) = arctan(z/z^) is the Gouy phase shift 



1 



■ arccos g 



(A4) 



which is determined by the 
-$mnp(0,0,- 



$ mnp (0,0,L/2) 



boundary condition 
L/2) = pn/2. The 
mode functions M^„ p (a;, y, z) are proportional to the elec- 
tric field and are normalized according to 



dx dy dzu° mnp {x,y,z)u m , n , pl (x,y,z) = 5 mm >5 nn ,S p 



Notice that the mode functions depends upon to, n, and 
p also implicitly, through the dependence of the wave 



vector k 



o/c. 



The thin membrane is a dielectric slab of thickness 
L d and index of refraction um, and when it is placed 
within the cavity the mode functions and their frequency 
change in a way which is dependent upon the position and 
orientation of the membrane with respect to the cavity. 
In order to minimize diffraction losses, the membrane 
should be placed orthogonal to the cavity axis and ex- 
actly at the cavity waist, so that the optical wavefront 
coincides as much as possible with the membrane sur- 
face. In the limiting case of a thin membrane, placed at 
Zo close to the waist, well within the Rayleigh range so 
that Zq ± L d /2 <C zr (see Fig. 1), the above Hermite- 
Gauss modes can be still used to describe the modified 
optical modes. In fact in this limit the new mode func- 
tions arc well approximated by 



(x, y, z) sin [$„ mp (x, y, z) + $ mrip (0, 0, L/2)} -L/2 < z < z - L d /2 

(x, y, z) sin [$ m „ p (x, y, z) + L d /2 < z < zo + L d /2 (A5) 

B mnp T mnp (x, y, z) sin [$ m „ p (x, y, z) - $ mrip (0, 0, L/2)] z + L d /2 < z < L/2, 



so that they automatically satisfies the null boundary 
condition at the mirrors. The dimensionlcss parameters 
Amnp, B mnp , C mnp , and ip mn p are determined by match- 
ing the solutions at the two membrane boundaries. The 
boundary conditions also determine the implicit equation 



for the new mode frequencies iu. 



mnp 



ch 



mnp : 



sin [k mnp L - 2{n + to + l)tp(L/2) + f3(k mnp )] 

)cos(2fc„ mp z ), (A6) 



where 



K(k) 



I) 5 



4n 2 M cot 2 (n M kL d ) + (n 2 M + l) 2 



(A7) 



is the membrane intensity reflection coefficient 45|, and 



kL cl 



arccos 



arcsm 



2n M cot(n M kL d ) 



\j^n 2 M cot 2 (n M kL d ) + (n 



M 



- 1 



— kL„ 



(A8) 



The membrane causes a frequency shift of each cavity 
mode which depends upon the position of the membrane 
along the cavity axis zq. However, this shift is always 
much smaller than the bare frequency, i.e., the solution 



of Eq. (|A6[) can be written as k Tl 



mnp 



with 
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\8k mnp \ < k^p, and using the fact that $o( fc m« P ) = 
2$ mnp (0, 0, L/2) = k^ np L - 2(n + m + 1)^(V 2 ) = 
one can derive the explicit expression of the frequency 



shift of a given cavity mode with unperturbed frequency 

, .0 _ _jL0 

^nmp nrap 



8u) r . 



1) 



(-lfy^ (fc° mp ) cos(2fc° r 



fi(k nrnp ) 



n Mk^ lmp L c i 



(A9) 



with [x] denoting the integer part of X. A selection of the 
modified cavity mode frequencies uj mnp = w^ np + 8uj nmp 
in the presence of the membrane are plotted versus the 
membrane position along the cavity axis zq in Fig. [Jj As 
predicted by Eq. (|A9[) they are periodic with a period 
equal to half wavelength. 

tem 00jP 

TEM 22 ,p-2 
TEM44„_ 4 
TEM n ,p., 
TEM 33i p_3 

. -TEM55,p-5 
-*TEM 12 ,p. 2 

TEMoi.p.1 
TEM 23iP . 3 

X 

TEMoo.p^ 

FIG. 7: (Color online) Frequency shift with respect to a 
driving laser with wavelength A = 1064 nm of some of the 
Hermite-Gauss modes versus the membrane position along 
the cavity axis zo. The plot refers to a cavity with length 
L = 9 cm, mirror radius of curvature R = 10 cm, and to a 
membrane with thickness Ld — 50 nm, and index of refraction 
tim — 2.0, which corresponds to a reflectivity \pfL ~ 0.4. 

The optical modes are therefore described by the "free 
optical" Hamiltonian 




H, 



opt 



I 



hujiajai, 



(A10) 



where uii 



are the frequencies of the optical modes 



in the presence of the membrane and I O {mnp} is a col- 
lective index. The creation operator aj creates a photon 
in the cavity mode with the mode function it mnp (x, y, z) 
of Eq. (|A5j) . and satisfies the boson commutation rule 



ai,a 



5w 



Considering the membrane motion means assuming 
that its position along the cavity axis oscillates in time, 
zq —> zq + z(x,y,t), where z(x,y,t) gives the membrane 
transverse deformation field. One typically assumes the 
high stress regime of a taut membrane, in which bending 
effects are negligible and one can use the classical wave 
equation to describe the membrane vibration. In such a 
limit the (normalized) vibrational normal modes of the 



square membrane with length D are 



<t>jk{x,y) 



D 



sin 



'(All) 

with j, k = 1, 2, . . ., \x\, \y\ < D/2. The vibrational mode 
with j, k indices oscillates in time with eigenfrequency 
Ojfc = (c s 7r / 'D)y/ j 2 + k 2 , where c s = y/T/a is the sound 
velocity, with T the surface tension and a the surface 
mass density. The membrane motion is then quantized 
by associating to each normal mode a bosonic annihila- 
tion operator hi (i O {jk} is the membrane collective 



index) such that 



bi,bl 



= 8, 



The free vibrational 



Hamiltonian can then be written as 



a?) 



(A12) 



where we have introduced the dimensionless position 
9« = m, + bi)/\/2 and momentum pi = (hi — b[)/"j2 
operators of each vibrational normal mode. The mem- 
brane deformation field operator can then be written in 
the equivalent forms 



z(x,y) = 



2aa 



(A13) 



^ * • ( 3' KX J 7T \ ■ ( kny kix 

— o, sin 1 sin 1 

2 J \ D 2 



D 



where m = <tD 2 /4 is the effective mass of each vibra- 
tional mode. 



1. The radiation pressure interaction Hamiltonian 

At lowest order in the membrane deformation field 
z(x, y), one can write the interaction Hamiltonian (in the 
Schrodinger picture) between the optical field in the cav- 
ity and the membrane motion as [46| 



Hi, 



D/2 



dx 



-D/2 



D/2 



-D/2 



dyP(x,y)z(x,y), (A14) 



where P(x, y) is the radiation pressure exerted by the 
cavity field on the membrane. A direct way of evaluat- 
ing this radiation pressure is to start from the Lorenz 
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force per unit volume f(x,y,z) = p(x,y, z)E(x,y, z) + 
J(x,y,z) x B(x 7 y,z) within the membrane [47j and in- 
tegrate it over the membrane length, that is 

P(x,y) = dzf z (x,y,z). (A15) 

Jz -L d /2 

The electric field E(x, y, z) is polarized parallel to the 
membrane and therefore its force contribution is balanced 
by the membrane support; moreover the current density 
J is due only to bound charges, J = dP/dt with P the 
polarization vector, and therefore the radiation pressure 
acting on the membrane can be written as 



P(x,y) =E Q (n 2 M -l) 



ZQ + - 



dz 



dE, 



z) x B(x,y,z) 
(A16)' 



The x axis can be taken along the direction of field po- 
larization and therefore the electric field operator within 
the cavity can be written (in the interaction picture) as 



E(x, y, z,t) = x 




ui(x,y,z)i 



aie~ iuit - aU iutt 



(A17) 

where ui(x, y, z) are the normalized mode functions spec- 
ified by Eqs. (|A5|) and we have used the collective index 
I <-> {mnp}. The membrane is placed close to the waist, 
well within the Rayleigh range zq ± L d /2 <C zr, and this 
allows to simplify considerably ui(x,y,z) and therefore 
E(x,y, z,t) within the membrane. In fact, one can take 
,w(z) ~ wq, ip(z) ~ 0, and R(z) ~ oo in Eq. (|A3|) , so 
that $i{x, y, z) ~ timMz, and by using Eq. (|A5p . one can 
rewrite Eq. (|A17[) within the membrane as 



E(x,y,z,t) = 



fkJi Cifi(x,y) . 



sm(n M kiz + tp t )i (aie 



where, from Eq. (|A2[I . 



fi(x,y) -.n^MHMwo] a 



w 



v / 7r 2"+ m - 1 n!m! 



x 2 + y 2 



(A18) 



(A19) 



is the normalized transverse pattern of the optical field within the membrane. The corresponding expression of the 
magnetic field within the membrane is given by 



St , ,n M -sr^ [f^Ji 
B(x,y,z) = y—^ ] / — 



haji CiTi(x,y) 



(A20) 



so that, inserting Eqs. (|A18|) and (|A20|) into Eq. (|A16|) and performing the integration, one gets the following general 
expression for the radiation pressure operator of the cavity field on the membrane (in the Schrodingcr picture) 



p(x, y ) = { <; 1)h 



yjujiw m fi{x, y)f m (x, y)h-i, m (aia m + aia^ + aja m + a\a\ n , 



(A21) 



with the dimensionlcss "longitudinal" factor A^. 



A 



kl ~\~ hji 



sin 



+ k m )L d 



sin[n M (fcz + k m )z + {(pi + <p m )] 



(A22) 



This factor depends upon the constants Ci and <pi associated with the cavity mode functions within the membrane 
(see Eq. ([A5[) ~). Notice that a similar term in A; >m with fc; + k m — > ki — k m gives no contribution when inserted into the 
sum of Eq. (|A21|) because it is antisymmetric with respect to I -<-> m exchanges. The cavity mode constants Ci and 
ipi depends upon the physical parameters of the cavity-membrane system, and their explicit expression is determined 
by the boundary conditions and by the normalization. After lengthy but straightforward calculations one gets 



A;, m = (n 2 M - l) A it 



ki + k n 

n(ki)n(k m 



n M {ki + k m )L d 



[1 -K{k m ) cos 2 {2k m zo)\ [I -TZ{h) cos 2 {2kizo)) 



1/4 



sin(2fc m z ) 
sm(n M kiL d ) 

s(fc,) l + y/l-S^m) 2 



sin[$ (fcQ-fc ; £ d ] 1 + ^1-sjh) 2 
sin [$ (&m) - k m L d ] l + ^1 - s(k m y 



1/2 



1 



(A23) 
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with 7Z(k) the reflectivity given by Eq. (|A7|) . and 



( — 



<Z> (k) = 2$ rmip (0, 0, L/2) = kL- 2(n + m + 1) arctan - — , (A24) 

. _ sm(2kz a )sm(n M kL d ) (aoh\ 
S[k > ~ sin [<M*0 ' ( } 



It is easy to verify from Eq. (|A23[) that when j = I the 
longitudinal factor Ajj becomes much simpler 

A u = M^izoj - 17( ^ h \ (oh v (A26) 

y 1 -^(fc;)cOS^(2fc;Z ) 

By comparing with Eq. (|A9I) . one can see that A;.; cx 
dui/dzQ, as it must be expected from the fact that the 
interaction Hamiltonian of Eq. (|A14[) is at first order in 
the membrane deformation field z(x,y,t). Eq. (|A26p is 
often useful in practice because one can safely approxi- 
mate Aj t i ~ Ay in the expression of the radiation pres- 
sure field of Eq. (|A21[) whenever | kj — ki\ <C h , which is 
verified when the driving laser is not too broad, so that 
only the cavity modes within a narrow frequency band 
are significantly populated. 

Eqs. (|A21I) - (|A25[) provide the general expression of 
the radiation pressure operator acting on the membrane 
modes which, inserted within Eq. (|A14[) yields the op- 
tomechanical interaction Hamiltonian at first order in 
the membrane deformation field z(x,y). Using the ex- 
plicit expression of Eq. (|A13|) . one finally gets the general 
expression of the radiation pressure optomechanical in- 
teraction for the membrane-in-the-middle configuration 




+a]a t + ajaj) , (A27) 

where we have defined the dimensionless transverse over- 
lap integral 6i,j,i between the j-th and l-th cavity modes 
and the i-th membrane vibrational mode, 

D ,D/2 ,D/2 

®i,o,i = ^" / dx / dy^i{x,y)fj(x,y)fi(x,y). 

Z J-D/2 J-D/2 

(A28) 

Due to the normalization of the involved mode functions, 
these overlaps satisfy |©i,j,/| < 1. 

Eq. (|A27p generalizes in various directions existing 
treatments of the radiation pressure interaction Hamilto- 
nian. In fact, it is possible to verify that it reproduces the 
radiation pressure interaction Hamiltonian of Ref. [23[ if 
one restricts to the longitudinal axis only and considers 
a one-dimensional cavity with a moving perfectly reflect- 
ing mirror TZ(k) = 1. It is also straightforward to verify 
that Eq. (|A27j) reproduces the expression of the radiation 
pressure operator of Ref. [43] in the free-space limit when 
there is no cavity. 



Since we are considering driven optical systems, the 
general interaction Hamiltonian of Eq. (|A27|) can be 
further simplified. In fact, at optical frequencies the 
counter-rotating terms ajcii + h.c have always negligible 
effects, and this holds true for the vacuum contribution 
associated with the Casimir force, as well As a con- 
sequence, for a typical membrane-in-the-middle system 
the multi-mode interaction Hamiltonian simplifies to 

Hint = ~ \ ~^ H ' ' /A •••'•'/'''>'/• ( A29 ) 

Therefore, the total Hamiltonian of the cavity-membrane 
system in the absence of driving is 

Htot = H opt + H V it + Hi n t (A30) 

with H opt given by Eq. (lATOl) . H vib given by Eq. (lA"T2l) . 
and by Eq. (|X29|> . 

Eq. (|A29[) shows that radiation pressure leads in gen- 
eral to a trilinear interaction describing all possible scat- 
tering processes in which a photon from one cavity mode 
is scattered to another one, mediated by the emission or 
absorption of a vibrational phonon. The optomechani- 
cal coupling rates associated with these scattering pro- 
cesses depends upon both cavity and membrane proper- 
ties, as well as the membrane position, as illustrated by 
Eq. |A"29|) together with Eqs. (jA"2"3|) and (|A"28|) . 

2. Reduction to the single-mode-single-resonator 
case 

One can restrict to a single cavity mode whenever the 
driving laser is tuned and aligned in order to be matched 
and to populate a given cavity mode only (here associ- 
ated with the annihilation operator a), and if scattering 
from the driven mode to other modes is negligible. Pho- 
ton scattering by the membrane can be neglected pro- 
vided that the membrane is not misaligned and not too 
far from the cavity waist, and if the motional sidebands 
of the driven mode do not overlap with nearby cavity 
modes. This latter condition is satisfied when the rele- 
vant mechanical frequencies fli are smaller than the typ- 
ical cavity mode separation, which is of the order of the 
free spectral range c/2L. Under these conditions, the 
interaction Hamiltonian of Eq. (|A29[) reduces to 

H int = -c^ay^giqi, (A31) 
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where gi is the vacuum optomechanical coupling rate be- 
tween the selected cavity mode and the i-th vibrational 
mode of the membrane, given by 

* - ^&- A »- < A32 » 

with ujq = cko the frequency of the driven cavity 
mode, 0^0,0 the transverse overlap of Eq. (|A28|) . an d 

A = sm(2k z ) {TZ(k )/ [1 - 1l(ko) cos 2 {2k z )] } 1/2 
[see Eq. (|A26|) ]. Eq. (|A31I) shows that the driven cav- 
ity mode is directly coupled to a collective membrane 
vibrational operator q e ff oc However, if the 

detection bandwidth is chosen so that it includes only 
a single, isolated, mechanical resonance with frequency 
f2 m , the collective coordinate q e ff becomes practically 
indistinguishable from the position operator q of the se- 
lected membrane mode, and a single mechanical res- 
onator description applies. When the term associated 
with the driving laser is included, the multimode Hamil- 
tonian of Eq. (|A30[) reduces just to single-mode-single- 
resonator Hamiltonian of Eq. ([1}, commonly adopted 
for the membrane-in-the-middle setup, provided that we 
identify the position-dependent frequency u>(q) with its 



linear expansion, 

u(q) =u) - gq. (A33) 

In fact, the present treatment is valid only when the 
membrane is placed out of the nodes and antinodes of 
the cavity field, where the linear expansion of Eq. (|A33p 
applies. If instead the membrane is placed exactly at a 
node or at an antinode, the resulting radiation pressure 
of the cavity field P(x, y) becomes equal to zero, imply- 
ing gi = 0. In such a case the optomechanical interaction 
becomes dispersive and is no more associated with radi- 
ation pressure. 



Appendix B: Stability conditions 

The stationary state of the cavity-membrane system 
under study is stable if and only if all the eigenvalues 
of the drift matrix A have negative real part. This con- 
dition can be re-expressed in equivalent form using the 
Routh-Hurwitz criteria [Hj], yielding the following three 
independent, nontrivial stability conditions: 



so 
8 1 



So 



fi m7m + 2K T (q s 



[irn + KT(q S )} 



n m GT > o, 

.2" 



2 lm n T {q s ) { \n T {q s ) 2 + (f2 OT - A) 2 ] \n T {q s ) 2 + (f2 m + A) 5 
+7m [[7m + 2K T (q s )} [n T {q s ) 2 + A 2 ] + 2K T (q s )n 2 n ] } + ACl m G 2 [ 7m + 2 KT (q s )Y 
+ {^Tm + 2kt(<Zs) [A 2 + [ 7m + K T {q s )} 2 ] } 

+n m (h lm - GT) [n m GT + 2n T {q s )tt m {h + O m ) + lm [K T (q s ) 2 + A 2 ]] 

) + 7™] 2 [h [K T (q s ) 2 + A 2 ] + K T (q s )GT] > 0, 
(fi™ + h) [Mis) 2 + A 2 ] - G 2 A + K T (q s )Gr > 0. 



(Bl) 



(B2) 
(B3) 



These stability conditions are more involved with respect 
to the standard optomechanical case without membrane 
absorption and mechanical frequency shift (h = T = 0), 
discussed in @, [3l|. In particular, there is the additional 
stability condition given by sq, which is instead always 
satisfied when h = T = 0: its violation corresponds to 
an "antistable" situation, i.e., when the vibrational reso- 
nance frequency becomes imaginary. We recall that the 
violation of the second condition, si < 0, corresponds to 



the emergence of a self-sustained oscillation regime where 
the effective mechanical damping rate vanishes. In this 
regime, the laser energy leaks into the sidebands and also 
feeds the coherent oscillations of the membrane vibra- 
tional mode. A complex multistable regime can emerge 
as described in [4!| . The violation of the second condi- 
tion S2 < instead corresponds to the emergence of a 
bistable behavior observed in (29j . 
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